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Exponential approximation based on the first order mean spherical approximation (FMSA) is applied to the 
study of the structure and thermodynamics of hard-core repulsive Yukawa fluids. The proposed theory utilizes 
an exponential enhancement of the analytical solution of the FMSA due to Tang and Lu [J. Chem. Phys., 1993, 99, 
9828] for the radial distribution function. From comparison with computer simulation data we have shown that 
at low density and low temperature conditions, where original FMSA theory fails, the FMSA-based exponential 
theory predicts a significant improvement. 
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1. Introduction 

Authors are pleased and honored to contribute to the festschrift dedicated to the sixtieth anniversary 
of Orest Pizio. Around thirty years ago one of us (AT) started his scientific career under Orest's super- 
vision. At that time we were dealing primarily with the electroneutral mixture of charged hard spheres 
that are widely known as a restricted primitive model of electrolytes. Such a model incorporates both 
Coulomb repulsion and Coulomb attraction. Repulsive interaction in the form of Yukawa potential, that 
is the subject of the present study, is quite frequently referred to as the screened Coulomb repulsion, 
and is widely used in modelling various real substances in statistical physics of soft matter. In particular, 
one can apply the Yukawa repulsion to the modelling of electrostatic interactions in one-component plas- 
mas |1], electrical double layer interactions in charge-stabilized colloidal dispersions and micellar 
surfactant solutions |:>]. etc. A brief but rather comprehensive review on the application of a repulsive 
Yukawa model in statistical physics can be found in the paper by Hopkins et al |6Q . Some of these applica- 
tions concern a pure repulsive Yukawa model, but very often the Yukawa repulsion is only a part of the 
potential functions used in modelling a wide class of the so-called complex fluids. Among others we wish 
to mention the core-softened fluid models, both homogeneous and those confined to a single surface, slit 
pore or even porous media. The field of core-softened fluid models is intensively developing nowadays 
and represents one of those in which Orest Pizio and his colleagues have recently taken interest 

Most generally, the repulsive Yukawa model is defined by the following potential function 

Y 

w(r) = -exp(-zr), (1) 
r 

where y > is the energy (coupling) parameter that characterizes the magnitude of repulsion while pa- 
rameter z controls the decay of repulsive interaction. This model has been extensively studied by means 
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of various methods of statistical physics, such as numerical solutions of Ornstein-Zernike integral equa- 
tion dEl EH, density functional theory (ljl[l2|] and computer simulations Blii HoHliHlil . In contrast 
to the attractive Yukawa potential, characterized by y < 0, the repulsive Yukawa model in general does 
not require the introduction of an extra core repulsion at short separations. However, looking for the 
analytical theory of the repulsive Yukawa fluid, attention was turned to the Waisman's solution of the 
mean-spherical approximation (MSA) for Yukawa interaction [17j. For reasons of the MSA solution tech- 
nique, the Yukawa potential is written in a form that slightly differs from that of equation (I), namely, 

f oo, r <a 

u{r) = \ eg (2) 
I — exp(-z[r-(T]), rscr, 

where g is the hard-core diameter, e = yexp(-zer) corresponds to the potential energy at a hard-core 
contact distance, r — g . Since the MSA solution is not sensitive to the sign of energy parameter e (or y ), 
it can be applicable to both repulsive and attractive Yukawa interactions. 

The MSA theory for Yukawa fluids being analytical in general, requires a solution of a set of rather 
complex algebraic equations. Although the original solution due to Waisman |17] was later on simplified 
by Blum and Hoye (3. Cummings and Smith |19]. Ginoza (2tJ and finally by Henderson et al {2lj and 
Scalise et al I22I1 . it still invokes a quadratic algebraic equation from which a unique physical root cor- 
responding to the coupling parameter (prototype of the Debye screening parameter) should be deduced. 
Recently, looking for further simplification of the MSA theory, Tang and Lu |23y have solved Ornstein- 
Zernike equation by means of perturbation theory and developed the so-called first-order mean spherical 
approximation (FMSA) theory. While the FMSA is a simpler theory, it is only slightly less accurate than 
the full MSA solution, and in contrast to the latter, is "solvable" at any thermodynamic conditions I24JI . Its 
simplicity, however, makes it also attractive for an employment in different related approaches, such as 
density functional theory HHIi^] or augmented perturbation theor y j27ll , etc. 

Similar to the MSA theory, the FMSA solution by Tang and Lu (2J is valid for Yukawa potential {2) 
in both e < and e > cases. However, in the present study our interest will lie in applying the FMSA 
theory to the case e > only, and we will refer to this model as the hard-core repulsive Yukawa (HCRY) 
fluid. It has been already shown |ipj,|28|] that FMSA theory is reasonably accurate for dense HCRY fluids 
and at high temperatures, i.e., when Yukawa repulsion is weak or dominated by a hard-core repulsion. 
However, the main difficulty of the MSA-like theories for repulsive potentials concerns the appearance of 
negative values of the radial distribution function at short separations in the case of the growing strength 
of repulsive interactions. Thus, although this has not been documented in literature, it is expected that 
FMSA theory will start to fail if temperature in the HCRY fluid becomes lower. In particular, in what 
follows we will show that such a deficiency of the FMSA theory becomes especially drastic at low temper- 
ature and low density conditions. It is obvious that in this region of density and temperature parameters 
the problems occur not only with the FMSA predictions for the radial distribution function but also for the 
thermodynamics of the HCRY fluid. A natural way of fixing such a flaw in any linear theory is to employ 
an exponential approximation for the radial distribution function l29l - [32ll . Besides ensuring positive val- 
ues of the radial distribution function at short separations, the exponential approximation also provides 
a correct asymptote of the radial distribution function at large distances. This significant improvement 
for the radial distribution function can be extended towards thermodynamics by means of the so-called 
energy route when the internal energy of the system is evaluated using its standard definition through 
the radial distribution function. The main purpose of the present study is to employ this strategy for the 
HCRY fluid based on the existing solution of the FMSA theory for Yukawa potential. 

In the following section 2 we will introduce and determine what we are referring to as the simplified 
exponential (SEXP) theory based on the FMSA solution due to Tang and Lu |23J]. This will be done for 
both the radial distribution function and for the thermodynamics. The results and their discussions are 
the subject of section 3. Since the low temperature and low density states have not been studied for the 
HCRY fluid so far, to verify the performance of the SEXP/FMSA theory we were forced to carry out the 
corresponding computer simulations. Thus, the necessary details of the canonical and constant pressure 
ensembles Monte Carlo techniques are briefly outlined as well. Our concluding remarks are summarized 
in the last section 4. To keep the paper simple and easy to understand, some auxiliary equations are 
moved to the appendix. 
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2. Exponential theory based on the FMSA solution 

The structure and thermodynamics of the HCRY fluid, defined according to equation (2}, will be stud- 
ied here within the exponential approximation based on the analytical solution of the first-order mean- 
spherical approximation (FMSA). The complete scheme and basic expressions of the FMSA solution for 
the Yukawa fluid are presented in separate papers by Tang and Lu (23 01- Here we recall only the equa- 
tions necessary for the radial distribution function. 



2.1. FMSA solution for the radial distribution function 

Like in any linear theory, the radial distribution function within the FMSA theory is represented as a 
sum of two terms, 

g™ SA {r) = g {r)-l3£g 1 {r), (3) 

where p = (^T)" 1 is the inverse of temperature, while go M is the radial distribution function of the 
reference system that is a fluid of pure hard spheres, and gi (r) is the distance dependent correction term. 
The FMSA theory due to Tang and Lu provides us with the Laplace transform G\{s) for the correction 
term OSSl 



Gi(s)= f rgi(r)e- sr dr = 



e - 1 

(4) 



where function Qo(s) refers to the Laplace transform of the Baxter factorization function resulting from 
the analytical solution of the Percus-Yevick (PY) approximation for pure hard spheres liill . 



S(s) + 12r?L(s)e" 

(1-77) 2 S 3 



with 77 = n pa 3 /6 being the packing fraction parameter, while L{s) and S(s) are polynomials. The equa- 
tion for the Laplace transform, Gq(s), of the radial distribution function of a fluid of pure hard spheres 
within the framework of the PY reads Ba, 13411 



ex 

/ 



L(s)e" 
(l-7?) 2 Qo(s) 



G {s)= rg (r)e- M dr = - - ? -. (6) 



For the purpose of the present study it is important to have rather compact expressions in order to 
invert the two functions, Gi (s) and Gq{s), into a real space. In the case of pure hard spheres, it reads (iH 

Hi 

^ 00 

go(r) = — — T [-l2T]) n Ca,n + l,n + l,rl(T-n-l) . (7) 
A similar inversion equation for the FMSA correction term is given by fiil liill 

1 (I-77) 4 00 



( 1 — 77) / V \ 

gitr)= £ {n + \){-\2r]) n D 6,n,n + 2,z<r, - - n- 1 . 



The functions C{n\,n2,n?,,rlo) and D{n\, 112,113, za,r/a) as well as the still undefined polynomials L(s) 
and S(s) are introduced in appendix. 



2.2. SEXP/FMSA approximation for the radial distribution function 

Equations (3)-GD represent the FMSA theory for the radial distribution function of the hard-core 
Yukawa fluid in general, i.e., for both the repulsive and attractive types of interaction. When evaluated 
numerically, the results predicted by equations {3)-l|8) are both qualitatively and quantitatively quite sim- 
ilar to those obtained within the complete MSA theory and are rather satisfactory in all practical applica- 
tions in the case of the attractive Yukawa interaction, i.e., when £ < 0. However, the same FMSA theory 
as well as the original MSA theory for the radial distribution function are not always that efficient in the 
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case of repulsive Yukawa interaction, i.e., for the HCRY fluid which is the subject of the present study. 
In particular, by lowering the temperature for the HCRY fluid, one is actually increasing the strength 
of repulsive interaction. For such conditions, the radial distribution function g(r) , being evaluated by 
means of equations (|3)-l|8), becomes negative at short distances r between particles. This is illustrated 
in figures [l] [2] and [3] where the results obtained within the FMSA theory are marked by dashed lines. 
Obviously, these physically incorrect results for the radial distribution function will also affect the FMSA 
predictions for the thermodynamics of the HCRY fluid. 




r/a r/a 

Figure 1 . (Color online) Radial distribution functions of the hard-core repulsive Yukawa fluid with za = 5 
at temperature T* = 0.125 and for several densities per 3 = 0.1; 0.2; 0.4 and 0.6 from the left-top to the 
right-bottom. Solid black lines denote the results of the FMSA theory, while the solid red lines mark the 
results of the SEXP/FMSA approximation. Symbols correspond to the MC simulation data. 

A natural way of fixing such a drawback of the linear theory is to employ an exponential (EXP) ap- 
proximation for the radial distribution function as it has been suggested within the MSA theory (illl . In 
the case of the FMSA theory, this approximation reads 

g SEXP (r) = go(r)exp[-/tegi(r)]. (9) 

Since FMSA theory itself is a simplified version of the more general MSA theory, the exponential approxi- 
mation based on the FMSA solution is usually referred to as the simplified exponential (SEXP) approxima- 
tion liill . By comparing equations (9) and {3} one can immediately conclude that g SEXP (r) is always posi- 
tive, which is due to the FMSA correction, /3cgi (r) , residing in the exponent. These results are shown by 
solid lines in figures[T][2]and[3] Additionally, exponential approximation l[9) results in the radial distribu- 
tion function that is asymptotically correct in the low density limit as well. Obviously, all these improved 
features of the exponential approximation remain valid within the original MSA theory. The reason for 
considering the SEXP/FMSA theory herein is due to the simplicity of the FMSA solution at nearly the same 
level of accuracy as the original MSA solution. This circumstance has caused a growing popularity of the 
FMSA theory itself. To illustrate this claim we are pointing out that in contrast to the MSA theory, the 
correction function gi (r) given by equation jg} does not depend on the temperature. Hence, there is very 
simple linear temperature dependence of the FMSA correction for the radial distribution function. This 
feature will be exploited herein below to write down the thermodynamics within the SEXP/FMSA theory. 
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2.3. Thermodynamics within the SEXP/FMSA approximation 

Thermodynamics within the SEXP/FMSA theory will be obtained herein through the so-called energy 
route 1301. Following this scheme, the derivation of thermodynamics starts with the evaluation of the 
internal energy E using its definition through the radial distribution function, 

oo 

pjj=2itppj drr 2 g s ™{r;P)u(r). (10) 
o 

This result is used to obtain the Helmholtz free energy in the form 

P oo p 

P^T = P^ + f d^'| = ^ + 12i7 / drr^r) f d/3'g SEXP faff) , (11) 
a o 

where Aq is the Helmholtz free energy of the fluid of hard spheres. Since functions go (r) and g\ (r) are 
independent of the temperature (see equation GD for g SEXP (r)), the inverse temperature integral in the 
second term of equation dill can be evaluated analytically, 

P 

I f d/3'g SEXP (r;/3') = go(r)exp[-/3egi(r)] V »o[-/3egi(r)], (12) 

where we introduced notation q>o (x) = [l -exp(-x)] I x. Then, it is more convenient to rewrite the re- 
sulting Helmholtz free energy expression in the form 

A- A id A FMSA A A 
P -=j8 +P , (13) 

N N N 

where we have introduced the conventional FMSA Helmholtz free energy contribution 

iFMSA A T C 

P— — =/3-^ + 1277j Arr 2 u{r) J d/3'g FMSA (r;/3') = a + a x + a 2 , (14) 

a 

with coefficients ao, a\ and a 2 for the HCRY fluid given by l^liill 

4n - 3ti 2 

ao = - 1 -y, (15) 



a 2 



oo P 

= 1277 [ drr 2 u{r) [ dff'go [r,p') = - 12 ^ L(z) (16) 
I { (l-r ? ) 2 Qo(z)z 2 

oo P 22 

= -\2t]eI Arr 2 u{r) f dp' gl {nP')p' = ^f-, (17) 



o 

and the correction due to the exponential approximation is 



oo 

/3^ = 1277/3 J drr 2 u(r)| g (r)exp[-/3egi(r)]<po[-/3£gi(r)]-go(r) + i/3£-gi(r)j. (18) 

a 

In contrast to equation (14) , the integral in equation < TT8t should be evaluated numerically by integrating 
along the radial distance, r, using the Gaussian quadratures. 

The pressure, P, can be computed using a standard thermodynamic relation, 

Pj = PV-fi±, (19) 
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where /i is the chemical potential. Within the SEXP/FMSA theory, the chemical potential can be repre- 
sented in the form, 

d ( A-A id \ d ( A ¥MSA \ d I AA\ fmsa 

This expression for chemical potential of the HCRY fluid is composed of two parts. The first part refers to 
the chemical potential within the FMSA theory, while the second one corresponds to the correction due 
to an exponential approximation for the radial distribution function. These two terms can be obtained 
from the following equations, 

jSm fmsa = f^o + ^i + ^] i + fl2> (21) 

\ dp dp dp j 

oo 

/3A/i = 1277/3 J"drr 2 u(r)j g (r)exp [-/tegi (r)] tp [~Pegi (r)] -go (r) + \p^gi w| 

OO 

+ 12P77/3 j drr 2 w(r) j -y^exp [-)3egi (r)] ?>o [-/teft (r)] 

£7 

« 3gi(r) r , r , 3g (r) 1 dgi(r)l 

+ fig ap go(r)exp[-jS£gi(r)Jyi[-j8egi(r)J- — + -/3e — j- , (22) 

where ip\ (x) = [l - x— exp (-x)] lx 2 , and five density derivatives were introduced, i.e., dag/dp, da\ldp, 
da,2ldp, dgo {r)ldp and dg\ {r)ldp . The first two of them could be straightforwardly evaluated using 
definitions (15)— (TT) for the quantities ao, a\ and «2» respectively. As for the next two, in accordance 
with equations and (8), one will be required to evaluate the density derivatives of the functions 
C{n\,ri2,i%,x) and D[ni,n2,ris,y,x). These derivatives are also straightforward and the details can be 
found in the appendix. 



3. Results and discussions 

We considered three models for the HCRY fluids with a fixed strength of repulsive energy, e = l, but 
three different ranges of repulsive interaction, i.e., characterized by three different values of the decay 
parameter, namely, za — 5, 3 and 1.8. The reduced temperature was fixed at T* - Tie = 0.125. While 
three values of the decay parameter z have been already used in related studies Q3, 13 of the HCRY 
fluids, such a low value of reduced temperature has not been explored so far in theoretical studies of this 
type of model fluids. At the same time, the reduced temperature T* = 0.125 and even lower temperatures 
are quite obvious in computer simulation studies of the phase diagrams of the HCRY fluid flEll- While 
considering such a low temperature, our main intention, on the one hand, is to clearly illustrate the 
problems experienced by linear theories in the description of purely repulsive fluids, and, on the other 
hand, to show an extent to which the application of exponential theories can improve the theoretical 
treatment of repulsive fluids. 

To test the theoretical predictions, we performed canonical ensemble and constant pressure ensemble 
Monte-Carlo (MC) simulations. To this end, the modified codes from the book of Frenkel and Smit (37I1 
were used and the radial distribution function g(r), the compressibility factor Z - pP/p and the internal 
energy /3E were calculated. In particular, the radial distribution function was obtained from canonical 
ensemble MC simulations for 256 particles in a cubic simulation box with periodic boundary conditions. 
All these simulations were performed for 500 equilibration steps and 1500 production steps. Each step 
consisted of 2500 attempts to displace a particle. On the other hand, constant pressure ensemble MC 
simulations were used to sample the density dependence of the pressure and internal energy. This set 
of simulations was performed for 256 particles with 2500 equilibration steps and 5000 production steps. 
Similarly, each step consisted of 10000 attempts to displace a particle and 100 attempts to change the box 
volume. 
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First, for each of three considered HCRY fluids we compare the radial distribution functions evalu- 
ated within the FMSA and SEXP/FMSA theories against MC simulations data. These results are shown in 
figures from[T]to[3]for the HCRY fluids characterized by za -5,3 and 1.8, respectively. Each part of these 
figures consists of a set of two theoretical (lines) and one computer simulation (symbols) results for one 
of four densities, namely, per 3 = 0.1, 0.2, 0.4 and 0.6. As it was already mentioned, the main problem of 
the radial distribution functions evaluated within the linear theories, such as MSA and FMSA, concerns 
an appearance of the non-physical negative values of g(r) at short distances r, when the strength of 
repulsion is increasing. Indeed, as we can see from figures [T| [2] and [3] this is the case at the considered 
reduced temperature T* = 0.125, where radial distribution functions shown by dashed lines (i.e., those 
that correspond to the FMSA theory) lie in the region of negative values. This region of negative values of 
g(r) is more notable and extends far beyond the contact distance rla = 1 for lower densities per 3 = 0.1 
and 0.2. We note that for these two densities, the distances r for which g(r) assume negative values, 
increase with an increase of the range of repulsion, i.e., going from za - 5, to za - 3 and 1.8. However, 
for density pa 3 - 0.4 and larger densities, the tendency changes to the opposite. Namely, the region of 
negative values of g(r) of the HCRY fluid characterized by za = 5 (see figure[l). starts to shrink going to 
fluids with za = 3 and 1.8. In the case of the highest density explored in this study, per 3 = 0.6, the radial 
distribution function g(r), evaluated from the FMSA theory, becomes even fully positive for za — 3 and 
1.8 (see figures[2]and[3), while it remains always negative for the case of za = 5. 



Table 1. Contact values of the radial distribution function of the HCRY fluids with a different range of 
repulsive interaction characterized by different values of the decay parameter z. Reduced temperature 
is fixed at T* =0.125. 



3 

per- 5 




za — 5 






za = 3 






zer = 1.8 




MC 


FMSA 


SEXP 


MC 


FMSA 


SEXP 


MC 


FMSA 


SEXP 


0.1 


0.003 


-5.943 


0.001 


0.001 


-5.549 


0.001 


0.007 


-5.154 


0.002 


0.2 


0.006 


-4.891 


0.003 


0.008 


-4.209 


0.005 


0.025 


-3.580 


0.010 


0.4 


0.028 


-2.794 


0.018 


0.081 


-1.817 


0.049 


0.237 


-1.057 


0.105 


0.6 


0.129 


-0.667 


0.108 


0.451 


0.307 


0.286 


0.914 


0.945 


0.541 


0.7 


0.030 


0.446 


0.241 


0.920 


1.340 


0.590 


1.521 


1.876 


1.009 


0.8 


0.687 


1.633 


0.511 


1.668 


2.414 


1.114 


2.449 


2.843 


1.712 



More information concerning the radial distribution functions can be found in table [T] where the 
results for contact values, g(r = a), are collected. Looking for the MC data one can see, that for each of 
three HCRY fluids the contact value g(r = a) increases when density increases. However, in the case of 
the HCRY fluid with za - 5 the radial distribution function at the contact always remains smaller than 
one, g(r = a) < 1, while in the case of two other HCRY fluids with za — 3 and 1.8 this is not the case. 
From the results presented in table[T]it follows that SEXP/FMSA theory results in g(r = a) values that 
are always slightly lower than similar MC data, while in the case of the FMSA theory this depends on the 
density: at low densities FMSA theory results in contact values g(r = a) that are lower than MC data, but 
at high densities it reverses, i.e., the MC data for g(r = a) are smaller than those from the FMSA theory. 
Thus, there is a crossover density which depends on the range of repulsion (i.e., depends on the value of 
parameter z ). The existence of this crossover density explains why for some HCRY fluid there is a range of 
densities, where the FMSA theory results in the contact values g(r - a) that are very close to the MC data. 
However, it is evident that on average in all considered cases the SEXP/FMSA approximation improves an 
agreement between theory and computer simulations, especially, at short distances. At intermediate and 
large distances, both FMSA and SEXP/FMSA theories are quite similar, and both slightly underestimate 
the ordering that takes place in the HCRY fluids. In particular, the oscillations of the radial distribution 
functions predicted by the FMSA and SEXP/FMSA theories show to be slightly weaker than those resulting 
from computer simulation data. 

Corresponding results for thermodynamic properties are presented in figures [4] [5] and [6] and col- 
lected in tables [2] [3] and H] In these figures and tables we show a comparison of the density dependence 
of the internal energy and compressibility factor for all three HCRY fluids, obtained from the FMSA and 
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Figure 4. (Color online) Internal energy (left plot), and compressibility factor (right plot) of the hard-core 
repulsive Yukawa fluid with za = 5 at temperature T* = 0. 125. Solid black lines denote the results of the 
FMSA theory, while the solid red lines mark the results of the SEXP approximation. Symbols correspond 
to the MC simulation data. 
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Figure 5. (Color online) The same as figure|4]but for za = 3. 
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Figure 6. (Color online) The same as figure|4]but for za = 1.8. 



23003-9 



S. Hlushak, A. Trokhymchuk 



Table 2. Compressibility factor Z - fiPl p and internal energy fiE of the HCRY fluid with za = 5 at re- 
duced temperature T* =0.125. 







pPIp 






pE 




3 

pa 


MC 


FMSA 


SEXP 


MC 


FMSA 


SEXP 


0.087 


1.840 


1.310 


1.872 


0.297 


-1.705 


0.333 


0.146 


2.748 


2.351 


2.799 


0.623 


-2.092 


0.721 


0.203 


3.948 


3.902 


3.992 


1.072 


-1.961 


1.246 


0.272 


5.892 


6.381 


5.835 


1.818 


-1.222 


2.082 


0.385 


10.377 


11.639 


9.853 


3.622 


1.132 


3.933 


0.494 


16.195 


17.744 


14.855 


6.056 


4.371 


6.182 


0.629 


25.444 


26.686 


22.939 


10.067 


9.331 


9.549 


0.866 


46.173 


47.493 


44.759 


19.139 


19.684 


17.191 


Table 3. Compressibility factor Z = 
reduced temperature T* =0.125. 


PPIp and internal energy /5E of the HCRY fluid with za 


= 3 at 






PPIp 






PE 




pa 3 


MC 


FMSA 


SEXP 


MC 


FMSA 


SEXP 


0.116 


3.450 


2.820 


3.338 


1.175 


-2.337 


1.128 


0.161 


4.971 


4.660 


4.757 


1.957 


-2.014 


1.861 


0.219 


7.321 


7.556 


6.938 


3.218 


-0.889 


3.037 


0.321 


12.448 


13.623 


11.688 


6.119 


2.538 


5.717 


0.427 


18.716 


20.439 


17.524 


9.836 


7.232 


9.074 


0.567 


28.207 


29.845 


26.642 


15.591 


14.233 


14.140 


0.817 


48.982 


49.399 


48.497 


27.545 


27.373 


24.965 


Table 4. Compressibility factor Z = 
reduced temperature T* =0.125. 


PPI p and internal energy fiE of the HCRY fluid with za = 


1.8 at 






PPIp 






PE 




pa 3 


MC 


FMSA 


SEXP 


MC 


FMSA 


SEXP 


0.055 


2.929 


1.887 


2.727 


1.136 


-2.277 


0.937 


0.088 


4.561 


3.455 


4.133 


2.176 


-2.495 


1.727 


0.122 


6.534 


5.685 


5.910 


3.486 


-1.991 


2.746 


0.169 


9.456 


9.244 


8.679 


5.517 


-0.419 


4.406 


0.258 


15.521 


16.759 


14.773 


9.951 


4.415 


8.331 


0.354 


22.614 


24.999 


22.096 


15.341 


11.095 


13.409 


0.486 


32.948 


35.792 


32.855 


23.341 


20.998 


21.155 


0.734 


54.504 


56.609 


56.117 


39.437 


39.313 


37.236 



SEXP/FMSA theories, against the corresponding MC simulation data. One can see that the above discussed 
failure of the FMSA theory for the HCRY fluids that results in the negative values of the radial distribu- 
tion function, is projected into the non-physical negative values of the internal energy for this type of 
fluids. Again we can see that SEXP approximation significantly improves theoretical predictions in all 
considered cases. This manifests itself in the internal energy for all three HCRY fluids being positive 
now, and, in general, more accurate than the one that follows from the FMSA theory. Quite similar, but 
less pronounced trend is observed for the compressibility factor. The SEXP/FMSA theory performs better 
everywhere, except at very high densities, where the internal energy obtained from the FMSA theory 
seems to be more accurate. However, this still remains questionable, at least for the case of HCRY fluid 
with za = 5, since according to the phase diagram for this fluid [ljj], the high densities (e.g., densities 
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pa 3 > 0.7) are in the vicinity of the transition region or even could be already associated with the solid 
phase. Similarly, the highest densities for all other HCRY fluids have to be treated with care. 



4. Conclusions 



Summarizing, a simplified exponential SEXP/FMSA theory was employed to describe the thermody- 
namics of the hard core repulsive Yukawa (HCRY) fluid. The theory is built up on an analytical solution 
of the FMSA theory for the radial distribution function due to Tang and Lu (2^1 . The main reason that has 
pursued us into this project were nonphysical negative values for the radial distribution function of the 
HCRY fluid that are predicted by the original FMSA theory at low temperatures. Obviously, this artifact of 
the FMSA theory will cause problems for the FMSA thermodynamics of the HCRY fluid as well. 

The results presented in this study are obtained for the low reduced temperature characterized by 
T* = 0.125. Such a low reduced temperature for the HCRY fluids has not been explored within the integral 
equation theory studies so far. However, there are computer simulation studies of the phase diagram of 
the HCRY fluid at such and even lower temperatures Q, 3- We have shown that the conventional FMSA 
theory at low density leads to the non-physical negative values of the radial distribution function at short 
distances and, consequently, to incorrect thermodynamics yielding the negative internal energy, but it 
seems to be reasonably accurate at higher densities. By contrast, the exponential approximation based 
on the FMSA solution significantly improves the theoretical treatment of the HCRY fluids. The SEXP/FMSA 
theory proves to be more successful for the HCRY fluids characterized by a more short-ranged repulsive 
interaction. In the case of long-ranged repulsive potentials, the exponential approximation performing 
reasonably well at low densities, tends to predict weaker oscillations of the radial distribution function 
at higher densities. 



A. Coefficients of the inversion of the radial distribution function and 
their density derivatives 



The other still undefined functions are as follows: 



C(ni,n 2 ,n 3 ,x) ■ 



E E 



(1+T7/2) " 2 
■B On, n 2 ,n 3 , i, a) + ZT^T S n 1 +n 2 ,3m 



D{m,n 2 ,n 3 ,y,x) 



2 m 

EE 

a=0/=l 



-yx _ t a x 



E 



J (-D ; '(?a + y) J V 



H{x), 



M J- 
(-iyB(ni,n2,n 3 ,i,a) (I + 77/2)" 2 



where 



B{n\,n 2 ,n 3 ,i,a) 



[ta + y] 1 



-1 m-im-i-ki 

E E 



e yX d„ 1+ n 2 ,3n 3 fH[x) , 



(1-7,) 2 " 3 

A{ni,n 2 ,ki,t a ) 



(« 3 - i + k 2 )\ (2« 3 - 1 - 1 - k x - k 2 )l , m-i-h 



and 



A;i!fc 2 !(«3-i-A;i-fc2)![(«3-l)!] 2 

a, 7. + \ v n 2 \(i + ni )\ ( 1 V, w-i jn+i-h 
A{m,n 2 ,ki,t a )= 2- T, ^TTT TTx 1+ o^ U + 27 7J f <* 

In these expressions H{x) is the Heaviside step function, while t a stands for the three (a = 0, 1,2) roots 
of the equation S (t) = with S (?) being the polynomial 

S [t) = (1 - j]f r 5 + 677 (1 - 77) t 2 + 1877 2 r- 1277(1 + 277) . 
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The roots are given by 



\-2 V + (2r 1 f) m (y + c a + y.c- a )] 



I-;? 



, a = 0,1,2, 



where 



f = 3 + 3tl-tf, y ± 



1+1 



2t7 4 



1/2 



1/3 



and c-e 



27ti/3 



Now the derivatives of the C and D functions read 



dC{ni,n2,n,3,x) 
dp 



E E 



df a fa;tD , . , t x dB(ni,n 2 ,n 3 ,i,a) 

x e " B{ni,n 2 ,n3,i,a) +e " 

dp dp 



7T 

+ - 

6 



1_\ «2 



2 (1-^) 



2n 3 



+ 2n 3 - 



2n 3 + l 



3n 3 f H{X), 



dD[ni,n 2 ,n 5 ,y,x) J ^ ^ 



dp 



U=oi~i[ <9p jtj ;! dp je y! 



(t a +y)' a=0i=l 



i-1 

E 



(-l)^(r a + y) J 'xJ 



;=o J 
dB {n\,ri2, « 3 , /dp .df a f? (m> n2,n^, i, a) 

[ta + y] 1 l ~3P [t a -y) i+l 

n 2 7ia 3 (1 + 77/2)" 2 - 1 n 3 ^ 3 (1+77/2) 



(-1) 



"2 



12 



,2n 3 



2^TT e yX5 «l + "2,3n3 



and 



dB[ni,ri2,n-i,i,a) n^Tia 3 
= — -: -B{ni,n2,n?,,i,a) 



dp 



3(1-7?)' 



1 tt3-i n-i-i-ki 

E E 

/c 1= fc 2 =0 



(1 _ j?)2 , i3 E E 



dA{n 1 ,n 2 ,k 1 ,t a ) /dp 
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dA{ni,n 2 ,ki,t a ) 

dp 



n 2 

E 



n 2 \{i + ni)\ 



{ri2 - i)na 3 ( 1 



1 + iJ (l + 2 J7 )' l2 - ( '- 1 C +I "- fa 
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In these equation 

dfq _ na^_ I 
dp ~ 6 I 



-77 I-77 



I 3»? J 3(277/)^ 1 M (57? d7? J 
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where 



dy± 
drj 



+ 



1 I8t] 3 4r] 4 df 



and 




6y2(l+27 7 4//2)H/ 2 P OT] 
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CnpoujeHe ei<cnoHeHTHe Ha6/ii/i>i<eHHfl fl/ia TepMOflMHaiviiKM 
TBepfloc<|)epHoro ruii/my 3 K)KaBiBCbi<MM BiflWTOBxyBaHHflM 

C. ryiyiuai^El a. Tpoxi/iMHyP 

tPaKy^bTeT xiMHHoT Ta 6iOMO/iei<y/iflpHoT iHXeHepiT, YHiBepcuTeT BaHflep6i/ibfla, HewBi/i, TeHHeci, CLUA 
iHCTHTyT cf>i3HKn KOHfleHCOBaHnx cHdeM, HAH yKpaiHM, 79011 ZlbBiB, yicpaiHa 

EiccnoHeHTHe Ha6/in>KeHHH Ha ocHOBi cepeflHbO-C(()epn4Horo Ha6/iwxeHHfl nepiuoro nopaflKy (FMSA) 3acTO- 

COByETbCfl flO BUBHeHHfl CTpyKTypH Ta TepMOflUHaMiKH TBepflOCCfiepHOrO nyiHHy 3 K)KaBiBCbKHM BiflLUTOBXy- 

BaHH^M. 3anponoHOBaHa Teopifl KopucTyeTbcn cfiaKTOM, mo icHye aHa/ii™4HwPi po3B'fl30K FMSA, a 3BH4awHe 
eKcnoHeHTHe Ha6/in>KeHHfl cyTTEBO Pioro noKpamye. LU/ihxom nopiBHHHHH 3 KOMn'KrrepHUM eKcnepuMeHTOM 
noKa3aHO, mo B 06/iacri Ma/inx rycTWH Ta HW3bKwx TeMnepaTyp, TaM fle Teopm cepeflHboro no/ia He npai^ioe, 
eKcnoHeHTHe Ha6/in>KeHHfl npwBOflHTb flo cyrreBoro BflOCKOHa/ieHHH TeopeTUHHUx nepefl6a4eHb. 

K/iHJMOBi c/ioBa: nni/iH 3 KiKaBiBCbKUM BifliuTOBxyBaHHRM, cepeflHbo-ccpeptmHe na6jw>KeHHR nepmoro 
nopnfiKy, eKcnoHeHTHe Ha6nw)KeHHfl, MoHre Kapjio MOftemoBaHHH 
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